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APPARATUS FOR SIMULATTNC THE DEFORMATION OF 
MATERIALS. NOTABLY OF SOFT POPY TISSUES 



lOQOl l This application claims Drioritv to Fre nch Patent Anolication 01Q1717. filed 
02/08/2001 including the soecification. drawings, claims and abstract, incorporated 
herein bv reference i n its entirety. 

BACKGROUND OF THE INVENTION 

10002 1 Tlie following description of the background of the invention is provided 
simply as an aid in understanding the invention and is not admitted to describe or 
constitute prior art to the invention. 

[00031 The invention relates to the simulation of deformation of materials, notably 
of soft body tissues. 

10004J^Usually, a display interface is provided for representing at least one object 

following a predetermined shape, displayed with a selected grid pattern^ A user 

interface provided with mechanical handling means allows the user to simulate one or 
more forces exerted globally on the object. 

1 00051 A memory is also provided for storing^mechanical parame^^^ 
of which the object in question consists^data as to the posit^^^^ 

at the vertices of a selected mesh, anc^^fqrce dato wh^^^^ 

position, stress to be exerted on the object. 

[0006 1 A computer ,fipg rati v e I y connec ted to the memory eyaluates new positions of 
the vertices, as a function of the stress exerted and the mechanical parameters of the 
material. In order to display the deformation and/or the return to a resting position of 
the moving object, the computer calculates the evolution of the positions of the 
vertices over time as a function of this stress. When it is used to simulate the 
deformation (or cutting) of body tissues, a qualified person (such as a surgeon) is able 
to follow the simulation of a deformation and/or a retum to a resting position. 
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J.000 7I Ln „tlie_ p^uJin shed application no. FR-97145 (I(Lbyjtbisj Vpplicanty the 

simulation of the defonnations of the object was based on a linear m odel, in the 

emb^dim.entjd.escn^^b 

has proved very promising in the fiel d of surgical assistance for simulating the cutting 
of soft b ody tissu es Hiver, skin, etc) . the simulation of larfi;e movements^ possibly 
non-linear, was not completel y natural. M_or_eoj^_er.ihls_a 

the reaction to a force exerted on lust anv type of deformable material, particularly in 
jJie-£asg-Q f large move njgjLt^ 

10008 1 Thus, simulated deformations with large shifts in position are not normally 
rel iable^TTie present inyention sets out to .improve the situati on . 



SUMMARY OF T¥ 



jVENTlON 



[0009| According to .one embodiment o f the invention, the computer comprises^a 

module capable of repeatedly determining, for each mesh element, the deviation 
between the current length of a side of the mesh and its resting length, and of storing 

the respectiye force data for each vertex of the mesh, and ^ module capab^^^ 

repeatedly determining, for each vertex, new positional data relating to this vertex as a 
function of the composition of the forces exerted thereon, data relating to at least one 
earlier position of the vertex and mechanical parameters of the material. 

400 1 01 Accordino to one embodimen t of the invention , the mesh pattern chosen 

comprises at least one triangular base. 

1001 1 J Thus, if the object is in the form of a hollow three-dimensional envelope, the 

mesh pattem chosen is preferably triangular, so that the computer is capable of 
determining the composition of the forces at each vertex of a triangle, as a function of 
the deviation between the current length of each side of the triangle and the length of 
this side at rest. 

10012] If the object has a solid three-dimensional shape, the mesh pattem chosen is 

preferably tetrahedral, so that the computer can determine the composition of the 
forces at each vertex of the tetrahedron, as a function of the deviation between the 
current length of each edge of the tetrahedron and the length of this edge at rest. 
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|0013 | A dvantageously, the memory is arranged so as to store, in conjunction with 
each mesh element, mechanical parameters of the material of the mesh, at least partly 
defined locally, particularly at the level of the mesh or elements thereof. 

J^001 4j According to another embodiment of the present invention, the computer is 

designed to determine said deviation between current length and resting length, jn 

order to estimate a derivative of the potential energy of deformation of each mesh, 
expressed as a function of a Green-St Venant tensor and mechanical coefficients 
inherent in the material in this mesh. 

I00151_The derivative of the potential energy thus expressed in relation to the 
position of the vertex makes it possible to obtain the force (or composition offerees) 
exerted thereon. 

10016 1 In estimating the derivatives of the potential energy at the vertices, thus 
expressed as a function of a Green-St Venant tensor, the computer is adapted to 
determine more precisely the differences between the squares of the current lengths 

and resting lengths of each side, jn order to detennine the composition of the forces 

mentioned above. 

10017! The mechanical parameters mentioned above advantageously comprise the 
Lame coefficients of the material in each mesh in question. 

1 001 81 The apparatus according to one embodiment of the invention, can thus be 
used to simulate large movements and/or non-linear elasticity, possibly of 
heterogeneous objects, formed of different materials. 

1 001 9 1 It is thus possible to simulate the deformation of a solid three-dimensional 
object, possibly enveloped in a hollow skin. In this case it will be advisable to apply a 
tetrahedral grid to the object and a triangular grid to the skin and to apply suitable 
mechanical parameters to each type of mesh. It is also possible to simulate the 
movements of soft objects (such as muscles) in the presence of rigid integrated 
objects (such as bones), as will be seen hereinafter. 

100201 The present invention may also take the form of a process for simulating the | 
deformation of objects, particularly of soft body tissues, comprising all or some of the 
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fol lowing step^^£££a) specifying a predeterm ined form for the obj ect at rest, ^b) 

making a grid image of said form, based on a selected pattern, preferably with a 

triangular base (triangle or tetrahedron), ^c^ storing the positi^^^ 

of the mesh in the object at rest, ;[d) simulating on^^ 

the object (or the movements of the object),Xej estjmatmg,^^^^ 

mechanical properties of the material and the positions of the preceding vertices, a 

resultant force on each vertex of each mesh, a^c^Q as a function 

on each vertex, determining the movement and, hence, the deformation and/or the 

return to a position of equilibrium of the vertices of the meshes. 

100211 Step £c) also involves aUributing indices to the elements of the meshes, 
within the scope of establishing a data structure, and, if appropriate, simplifying the 
structure of the object, as will be seen hereinafter. Steps (e) and £f) are based on a 
model developed by the Applicant which will be described in detail hereinafter. 

10022L ]t is to be understood that both the foregoi_nR_^cne ral description and the 
fbllQwinKjlelajJejLdes^ and explanatory onlv> and are not 

restrictive of the inv ention as claimed. 

BRIEF DESCRIPTION OF THE DRAWINGS 

j;0023J Further features and advantages of the invention will become apparent from 

studying the detailed description that follows and the accompanying drawing^ 

described below. 

|0024| Figure I diagrammaticallv shows an apparatus according to one embodiment 
of the inventioHe 

40025 1 Figures 2 A and 2B diagrammatical ly show triangular meshes, in the initial 
and deformed state, respectively, with notations for the points of the triangle 
subsequently used in estimating the composition of the forces at each vertex of the 
trianglq^ 

40026 1 Figure 3 A is a diagram of a tetrahedron which helps in the understanding of 
the notations used for a point PJ of the tetrahedron and the normal vector nj of the 
opposite foce^ 
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J0027| Figures 3B to 3D are diagrams of tetrahedrons which, together with the 

abovementioned £igu^ 3 A, help to show how the sums on the edges and pairs of 

edges of a tetrahedron are calculated, in order to evaluate a force to which a vertex Pp 
of the tetrahedron is subjected^ 

J 0028 1 Fis sures 4 A and 4B are djagrams to assist in the understanding of the 

attribution of indices in different triangular meshe^ 



1002 91 Fj«iires 5 A and 5B are diagrams of triangular mo^^^ comprising, 

respectively, a current vertex (represented by a square) and a virtual vertex 
(represented by a disc), before simplification, and the same gridded shape after 
simplificatioi\. .._ 

^003 01 Figure 6 diagrammatica^^^^ an example of a noncon 

mesl\. .__ 

4.003 H Figure 7 shows a diagrammatic flow chart of the pr9@"ess of the calculations 
of forces and positions of the vertices representing the object. 

DETAILED DE S-CRIZPm 

100321 The present description will refer to formulae and equations which appear in 
the annex. 

10033] The drawings, and the aiinex, essentially contain elements of a definite 

nature. They can therefore not only be used to gain a fuller understanding of the 
description, but may also contribute to the defmition of the invention, where 
appropriate. 

10034] Reference is made first of all to Figure 1 in which the apparatus is shown in 

the embodiment described in the form of a computer (PC). The computer comprises a 
permanent memory MEM (for example of the ROM type) in which are stored 
pro grams for calculating the positi^^^^ and the resultant forces at 

these vertices, as well as a processor pperatiyel y connected to a working memory 
MT (for example of the RAM type) for carrying out these calculations. 
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10035 1 T he abovementioned modules of the computer described hereinbefore may, 
if necessary, be in the form of gograms stored in the memory MEM, or may be 
loaded into the working memory MT, from reading a data carrier in an external 
peripheral (not shown) such as a disc drive or a CD-ROM reader, or received via a 
communications link, e.g. a MODEM. Of course, such 4 3iQgymis j^ppgar jo be 
particularly important means for carrying out the invention. Thus, the present 
invention also relates to jaroarams of this kind, possibly on a data carrier of the type 
mentioned above (disc, CD-ROM, etc), detailed embodiments of which are described 
hereinafter. 

100361 The computer further comprises a monitor VISU connected to the processor 
[iP via a graphic interface IG, the unit consisting of VISU and IG (referred to as the 
"display interface" hereinbefore) providing a display of an object the deformation of 
which is being simulated. The computer further comprises a user interface RJ 
connected to mechanical handling means such as a keyboard CLA or a mouse MO, 
enabling a user to simulate a force exerted globally on the object, or a change of 
position of this object, with initial deformation. 

X0037] ^The shape of the object to be displayed in grid form may be a surface shown 

in three dimensions. In this case, the mesh is triangular (figure 2A). The shape of the 
object to be displayed in grid form may be a three-dimensional solid and in this case 
the mesh is tetrahedral figure 231 If appropriate, a grid module is stored in the 
memory MEM and can be implemented by the processor joP, co-operating with the 
graphic interface IG, to display a shape of the grid- covered object on the monitor 
VISU. 

10038] ^The following is a description of a model, for example purposes only. 

developed by the Applicant and applied to an apparatus for simulating deformation. 
This model uses a resting form without mathematical induction in relation to a resting 
position (mathematical induction by shifts which is normally the basis of conventional 
models). 

J^0039J Let us consider first of all a surface shown in three dimensions, gridded with 

a triangular mesh pattern. A triangle TO = (PO, PI, P2) of the resting form, with edges 
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ihJhh)' This triangle (on the left in figure 2A) undergoes a transformation O to 
give the triangle T = (QO, Ql, Q2) (on the right in figure 2B). 

100401 The potential elastic energy of this deformation is expressed here from the 
Cauchy-Green C deformation tensor (formula Fl of the Annex). Looking at the 
finished element obtained by associating a set of linear interpolation functions with 
this triangle, the deformation gradient is expressed at every point X on the triangle by 
the formula F2, wherein A is the area of the triangle TO, nj is the unit vector 
perpendicular to the segment (Pj+1, Pj+2), and the aj *s are shape vectors of the 
triangle. 

10041j The deformation tensor C is then given by the expression of formula F3. 

Advantageously, this deformation tensor is invariant by translation and therefore does 
not depend on the choice of the origin of the reference point. Preferably, the origin of 
the reference point is placed in the centre of the circle circumscribed around the 
deformed triangle T, of radius r (figure 2A) , which allows the tensor C to be 
expressed in the form given in formula F4, where Lp is the length of the deformed 
edge joining the vertices Qj and Qk. 

1 0042 1 T he expression of the potential deformation energy (given by formula F5 in 
the annex) uses those of the trace of the tensor C (given by formula F6) and of the 
trace of the tensor C2 (which will be described hereinafter). 

10043J Starting from the classic equations in the triangle TO (cosine rule), which are 

expressed by the formula F7, the scalar products between the directing vectors are 
obtained, given by equation Rl in the annex. The respective traces of the tensors C 
and E are thus written according to formulae 8 and 9 in the annex. This term depends 
only on the variations in length of the edges of the triangle and its resting 
shape(represented by the magnitudes Ki). 

[00441 From the expression of tr (formula FIO), trE^ is deduced and, finally, the 
potential elastic energy of deformation of the triangle is given by equation R2 in the 
annex. 

100451 The elastic force applied by the triangle at one of its vertices is obtained by 
deriving the elastic energy relative to the position of this vertex. In the embodiment 
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shown in Figure 2A, the vertex PO has QO for its deformed position. Deriving the 
equation R2 in relation to QO, the first step is preferably to derive the terms L^p in 
order to obtain the expression of the force according to formula Fll, which is the 
equation of the elastic force exerted by the triangle at each of its vertices. This force 
is expressed solely as a function of the variations in length of the sides of the triangle 
(Li^ - li^) , the mechanical parameters of the material (A. and \x) and the shape of the 
triangle at rest (given by the terms Ki). 

100461 Advantageously, this force implicitly contains a property of rotationaly 
invariance as all its terms remain zero as long as the lengths of the sides of the 
triangle do not change, i.e. as long as the triangle is not deformed. Thus, the 
calculations of movements can be simplified by virtue of the fact that it is no longer 
necessary to compensate for rigid transformations, as invariance by rigid 
transformation is intrinsic to the formulation of the deformable model. 

100471 Extending the foregoing to a solid three-dimensional shape (having a 
volume), gridded with a mesh of tetrahedrons, by formulating the elasticity according 
to Saint Venant-Kirchhoff it is possible to express the gradient of deformation 
(formula F12), as well as the Cauchy-Green deformation tensor (formula F13) , where 
Aj is the area of the face opposite the vertex Pj, nj is the perpendicular to this face and 
V is the volume of the tetrahedron (Figure 3 A). 

10048} By placing the centre of the reference point at the centre of the sphere 
circumscribed around the tetrahedron T, it is possible to express the scalar product of 
the position vectors of two vertices of the tetrahedron as a function of the length of the 
edge that joins them and the radius r of the sphere, thus making it possible to write the 
deformation tensor according to formula F14. 

IQ049I The elastic deformation energy is expressed as a function of the Green-Saint 
Venant deformation tensor (bearing the reference E), but in a 3D grid (formula F15). 
The expression of the trace of the tensor C is also expressed with the scalar products 
(aj.ak), using the cosine rule for dihedral angles and as a function of the lengths of 
the edges of the tetrahedron (formulae 16 and 17). Consequently, the trace of is 
written according to formula PI 8. Finally, the Applicant has shown that the elastic 
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deformation energy is expressed by the formula F19 of the annex (complemented by 
the formula F20). 

100501_The expression of the elastic force exerted by the tetrahedron on one of its 
vertices is obtained by deriving the elastic energy relative to the position of this 
vertex. A formulation which is easy to use is expressed by formula F21 of the annexe 
and is interpreted using Figures 3A to 3D. 

100511_In formula F2I 

in the second line, k and m are the complementary indices of j and p, 

in the last line, m is the complementary index of j, k and p, and : Kjkmp = 

KjpKkm KjmKkp'KjkKnift^ 

100521 The total on the edges (j\p) indicates that it is necessary to travel over the 
three edges adjacent to the vertex Pp (on w^hich the force is exerted) whereas the total 
on the edges (j,k) (indices different from p) involves the three edges located on the 
face opposite the vertex Pp. 

Referring to Figures 3 A to 3D: 

- the total on the edges (j,p) involves three terms (Figure 3 A), 

- the total on the pairs of edges ((j,p), (k,p)) involves three terms (Figure 3B), 

- the total on the pairs of edges (0,p), 0*^)) involves six terms (Figure 3C) and 

- the total on the pairs of edges (Q,k), (m,p)) involves three terms (Figure 3D). 

[00531 Once again, there is rotational invariance of the above expressions, for a 
solid three-dimensional form, which advantageously makes it possible to avoid the 
appearance of unrealistic distortions when part of the object is rotated too far. This 
invariance also extends to movements of translation. Thus, the object may undergo 
any overall movement of translation and rotation while retaining zero internal energy 
and stresses. The Applicant has verified that in this case the rigid component of the 
deformation is zero. 

10054 1 M oreover, the discretization of the object using finished elements comes to 
approximate the elastic energy over a set of sub-domains defined by simple geometric 
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shapes. The above model thus suggests, for each sub-domain, a formulation of the 
energy which only depends on the change in shape of this domain. 

100551 Advantageously, displaying a simulated deformation while considering an 
initial resting shape (rather than initial resting positions) appears to be more natural 
than the deformations simulated with known models. Moreover, considering an 
initial resting shape avoids the problems of accuracy that arise when the model moves 
substantially away from its resting position. 

1005 6 J The following is a description of the preliminaries to using the 

abovementioned computer modules to determine the forces at each vertex of the 
chosen mesh. 

f 0057 1 First of all, a data structure should be defined. The Applicant has first dealt 
with the difficulty of attributing suitable indices to the elements of each mesh. 

10058j For triangular patterns, the data structure corresponding to triangulation 

consists of a list of triangles, a list of sides and a list of vertices: 

- for each vertex, an adjacent side is determined, 

- for each side, a pointer is determined on the two adjacent triangles, and the two 
vertices which form the side are determined, 

- for each triangle, the pointers are determined on the three vertices that form the 
triangle, the three sides of the triangle are determined, the three triangles adjacent to 
each side are determined, and three values giving the index of adjacency of the 
current triangle for each adjacent triangle are determined. 

i00S9J The orientation of the three vertices in the triangle is carried out as follows. 

X0060J For the grid £is a whole to be correctly oriented, each edge or side must be 

traveled over in two different directions when each of its two adjacent triangles is 
traveled over. Referring to Figure 4A, the side Ai of index i is the side opposite the 
vertex Si of index i in this triangle. The adjacent triangle to index i is the /ieiohborin^ 
triangle TV, adjacent to the side Ai. Finally, each triangle always has three adjacent 
triangles. In order to represent surfaces with borders (a plane has one border, a 
cylinder has two borders) the notion of so-called "virtual" vertices and triangles is 
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used (Figure 4B) . The virtual triangles are not displayed and the position of the 
virtual vertex is of no consequence to the geometry of the grid. 

100611 An example of index attribution will now be given by way of example, in 
C-H- language. 

i0062]^If the adjacent triangle number i is called neighbor[i], then the adjacency 
index "^whichNeighborfi]" is defined by A'^ J.^struction such as 

Neighbor[i] — ^>neighbor[whichNeighbor[i]]=this ; 

10063 J For reasons connected with a tetrahedral grid, the data structure of a 

tetrahedral volume consists of a list of tetrahedrons, a list of triangles, a list of edges 
and a list of vertices, 

J0064 i For each vertex, its position and an adjacent tetrahedron are stored. 

I006S J For each edge, the pointers on the two adjacent vertices and a pointer on a 

triangle adjacent to this edge are stored. 

100661 For each triangle, the pointers on the three vertices that form the triangle and 
the pointers on the two tetrahedrons adjacent to the triangle are stored. 

100671 For each tetrahedron: 

- the pointers on the four vertices that form the tetrahedron, 

- the pointers on the six edges of the tetrahedron, 

- the pointers on the four triangles of the tetrahedron, 
and 

-the pointers on the four adjacent tetrahedrons opposite each of the four vertices are 
stored. 

1 0068 1 E ach tetrahedron is oriented. Starting fi-om the orientation of the four 
vertices ABCD, the operator seeks to orient the four triangles : BCD, ADC, ABD, 
ACB. To ensure correct orientation of the grid, each triangle has to be traveled over 
in two different directions by traveling over each of its two adjacent tetrahedrons. 
Finally, each tetrahedron is always adjacent to four tetrahedrons. The notion of so- 
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called "virtual" vertices and "virtual tetrahedrons" is introduced to take account of the 
borders of the grid. 

i0069J„In order to simplify the structure of the deformed object, the grid is cut up by 
removing material (more precisely a vertex Sr, although this vertex cannot be a virtual 
vertex). 

100701 In the case of a triangular mesh, if this vertex is not adjacent to a virtual 
vertex, this vertex is regarded as being virtual like all these adjacent triangles. If the 
vertex is already adjacent to a virtual vertex, the side which connects this point to the 
virtual vertex is removed, and all the triangles adjacent to the virtual vertex are forced 
to be virtual (Figures 5A and 5B), 

100711 With a tetrahedral grid pattern, the procedure involves removing a set of 
associated tetrahedrons. Once the tetreihedrons have been removed the grid must 
remain true, as described hereinafter. For this, it may be advantageous to remove 
more tetrahedrons than those present in the set of tetrahedrons mentioned at the 
beginning. Virtual vertices may also be added or removed when creating associated 
components or removing the borders. 

100721 The Applicant has confronted the problem of ensuring conformity of the 
grid, a definition of which is provided hereinafter. 

1 0073 ! F or a triangular mesh to be true, each side must belong to more than two 

triangles and the topological nei ghborhood of each vertex must remain 

homeomorphous with a circle (or with a semicircle), i.e. if two triangles share the 
same vertex, there must be a set of triangles adjacent to each other which make it 
possible to go from one triangle to another. Referring to Figure 6, the triangular mesh 
shown is not true, as the two triangles have an intersection along a vertex, while there 
are no adjacent triangles making it possible to go from one of these two triangles to 
the other. 

1 00741 For a tetrahedral grid to be true, each triangle must belong to more than two 

tetrahedrons and the topological /lei ghborhood of each vertex must be 

homeomorphous with a sphere (or a hemisphere if it is on the border). Therefore, if 
two tetrahedrons share the same vertex, it is essential that there is a set of tetrahedrons 
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adjacent to one another which make it possible to pass from one tetrahedron to 
another. 

10075 ] The equation that governs the movement of the linear elastic model is the 

Newtonian differential equation EI in the annex. This equation is connected with the 
differential equations encountered in the mechanics of continuous media (equation 
E2). 

£0076 1 According to the theory of the finished elements, the mass matrix M and the 
attenuation matrix C are hollow and depend on the properties of each tetrahedron. 
Moreover, these matrices M and C may be considered here as being diagonal, i.e. the 
mass and the effects of attenuation are concentrated on the vertices of the mesh. This 
simpliflcation, known as mass-lumping, uncouples the movements of each of the 
nodes and thus makes it possible to write the previous equation E2 as a set of 
independent differential equations of the type El, with an equation for each vertex. 

100771 The choice of an explicit integration diagram leads to the expression of the 
elastic forces at the time t, to calculate the new position of the vertex at the time t+1, 
with the equation E3 in the annex. Advantageously, this equation E3 applies to all 
types of force, both linear and otherwise. 

10078J The explicit diagrams are only conditionally stable and therefore converge 

less rapidly than the implicit diagrams. 

[00791 Another difficulty arises : for the diagram to be stable, the time step used for 
the digital integration should be sufficiently short. The critical time step depends on 
the greatest intrinsic value of the rigidity matrix and of the local values of the mass 
and of the attenuation. To optimize the time step, the mass matrix is fixed to be equal 
to identity, and the optimization of the local values of attenuation, so as to limit the 
risks of oscillations, have the effect of improving the packaging of the linear system, 
and hence of establishing its resolution. 

i00801_Reference will now be made to Figure 7, wherein a flow chart 
diagrammatically shows the different treatment steps applied to a triangulated three- 
dimensional surface, by way of example. 
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lOOSlj The treatment begins with a first step 10, in which the structure of the data 

described above is applied, with a suitable attribution of indices. Then in step 12 
consideration is given to each triangle Tn, for which the Lam6 coefficients Xn and j^n 
are sought, as well as the area of the triangle A (step 14). For each side of the triangle 
(j being an integer from 0 to 2), the difference between the squares of the lengths of 
the sides is calculated, between the current length Ij and the previous length Lj of the 
side(step 16). From these differences in the lengths of the sides squared, the forces 
exerted at each vertex Sj of the triangle Tn are deduced (force Fj" calculated in step 18 
of the flow chart). Moving on to the next triangles, the forces are integrated over all 
the triangles by looping on the index n (step 20). 

100821 In steps 30 and 32, the new positions Q'n of the vertices Sj of the deformed 
triangles are determined, from their previous position Qn and the forces F". By 
looping on to the index n in step 36, integration is carried out over all the deformed 
objects. The new positions of the vertices Sj are preferably estimated by using the 
equations El, E2 and E3 mentioned hereinbefore. 

10083J In the flow chart in Figure 7, the two respective integrations (of the forces Fj" 

and of the positions of the vertices Sj in the deformed objects) are shown separately. 
Of course, it may be that some of the steps 12 to 20 may be carried out at the same 
time as or after some of the steps 30 to 36. This choice depends on the structure of 
the two modules described hereinbefore, contained in the abovementioned computer, 
and, if applicable, on the data structure. 

10084 ! A test 42 is preferably provided to find out whether the object has fully 
returned to its resting position, i.e. : is there an edge in the object the length of which 
differs substantially from the previous one? 

[00851 Ifthe computer does not find such an edge as it travels over the triangles, the 
object has returned to its initial resting shape and the treatment is stopped at 40. 

i0086J^Otherwise, an image of the object is displayed by the monitor VISU 
mentioned above and a time loop (step 44) is incremented, to pass to the next time 
and show a progressive return to a resting shape of the object, using equation E3, if 
appropriate. This time loop can be initiated after the first starting step 10. 
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[00871 Of course, the invention is not limited to the embodiment described 
previously by way of example, but also extends to other alternative embodiments 
defmed within the scope of the claims that follow. 

|0088) Thus, it will be understood that the invention is not limited to an application 
to simulating the deformation of soft body tissues. As described previously, the shift 
in position of a rigid object (for example a bone attached to a muscle) integral with 
the soft tissue can also be simulated. In some cases it may be appropriate to use a 
tensor (of inherently known expression) rather than a deformation tensor of the type 
described above. 

10089J ^The formula for the deformation energy W, in both a triangular and a 

tetrahedral grid, is provided hereinbefore by way of example and depends on the 
material of which the object whose deformation is being simulated consists. 
However, it is generally expressed by a combination of the traces of the tensor E and 
the tensor E^. In these particular cases, only the trace of the tensor E (or the tensor E^) 
may be used. 

1009 01 The foregoing description of a preferred embodiment of the invention has 
been presented for purposes of illustration and description. It is not intended to be 
exhaustivejo^rtoli mit the i nvention to the precise form dls_closed> ancLmodificatiq ns 
and variations are possible in light of the above teaching or may be acquired from 
practice of the invention. The embodiment was chosen andj3escjibejlin.^fe 
eAI3Ldnjtke_pjindp^^^^^^ of the invention and as a p ractical application to enable one 
skilled in the art to utilize the inven tion in various embodiments and with various 
modification are suited to the parti cular use contem plated . It is intended that the 
scope of the invention be defined bv the claims appended hereto and their equivalents. 
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